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OO Introduction 

^ jThe following fact in elementary complex analysis is well known. 

■ 

■ Theorem A. Every harmonic function u: {\z\ < 1} — t- C 0/ a complex variable is 
expressible as the sum of a holomorphic function and an antiholomorphic function. 

\ This principle has many uses. In particular, when a holomorphic solution to a problem 
l^is sought and a first attempt is made as a harmonic function, one may "throw away" the 
''^antiholomorphic part to obtain a holomorphic approximation. A classical example of this 
. .^principle in conformal mapping theory is found in the method of Fornberg [10] , in which a 
guess as to the boundary values of the mapping is expressed on the boundary of the unit disk 
as a Fourier series, whose coefficients give the sum of a power series in z (positive powers 
of e*^) and in z (negative powers) in the interior of the disk. The sum of these two is a 
harmonic function, whose antiholomorphic part is discarded in the algorithm. 

Theorem A holds in many generalizations of the field C of complex numbers, for example 
monogenic (hyperholomorphic) functions on quaternions [28] or on Clifford algebras [H [T3]. 
It also holds for monogenic functions from to EI [5l [12] . 

In this paper we show (Theorem 13.11) that the natural generalization of Theorem A does 
not hold for monogenic functions from M.^ to M.^. Therefore, the class of harmonic functions 
which have no such decomposition is new. We will use the term "contragenic" for harmonic 
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functions which are orthogonal to the monogenic and the antimonogenic functions in the 
sense of L2. We beheve they are of interest because of the relevance of monogenic functions 
in to physical systems pi]. 

Thus we initiate here the study of contragenic functions. The precise definitions for 
monogenic functions and related notions in are specified in section [H based on work of 
[2l O UHl [17]. In section [2] we summarize the necessary facts about the standard basis for 
homogeneous monogenic polynomials in the unit ball of and related spaces. We also 
calculate an orthonormal basis for the vector parts of the monogenic functions. In section [3] 
we prove the existence of contragenic functions and prove some of their basic properties. In 
particular we derive a Bergman kernel which annihilates precisely the contragenics. Finally, 
in section H] we give an explicit construction for a graded basis for the space of contragenic 
functions. 

1 Monogenic functions in 
1.1 Notation 

We will use fairly standard notation for the skew field of quaternions HI = = {x = 
XqCq + XiCi + ^262 + 3:363, Xj G M}. Here Cq = 1 is the unit and the multiphcation is 
determined by e^^ = —1 and Cj^Cj^ = ±6^3 where {ji, j2, js} is the set of indices {1, 2, 3} and 
the + sign is taken precisely when the cyclic order matches that of 1,2,3. See texts such 
as [T2I [T5l [28] for further general information. The M-subspace of quaternions x G H such 
that 63 = will be denoted © {0} or more simply just when there is no danger of 
confusion. 

Let us write dj = d/dxj for j = 0, 1, 2, 3. There is a great deal of literature (see [28] and 
references in ^) on the Cauchy-Riemann type differential operators operators 



which act both from the left and from the right on differentiable functions / = /o + /iCi + 
12^2 + fs^s defined in open subsets of H. The operator is a generalization of the operator 
d/dz on which complex analysis is based: functions for which -De/ = (resp. fD^ = 0) 
are variously called left (resp. right) Fueter-regular, monogenic, hyperholomorphic, among 
others; occasionally the roles of and Dm are interchanged in the terminology. 

In recent years some work has been done |2l [12] on the analogous functions from M.^ to 
and from to M.^, with a view to expressing and studying operators relevant to physics. 
Relatively little has been done for functions from to M^; in this regard we mention [9[ [TS]. 
This setting is particularly interesting because even though such functions are not conserved 
under multiplication by elements of © {0}, i.e., the algebraic structure of an algebra or 
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ring is lost, in many ways the monogenic functions behave more hke standard holomorphic 
functions in C (cf. Proposition II. ip . To be precise, for x = Xq + Xiei + a;2e2 G M^, let us write 
Sex = Xq, Vecx = X = XiCi + 0:262, x = x — x. Let f2 C be an open set, and T-Lr{Q) the 
space of real-valued harmonic functions defined in fl. We consider the (Moisil-Theodorescu 
type) operators D and D, defined by 

d = 8161 + 8262, D = do-d, D = do + d, (1.1) 

and define the set of (left-)monogenic (or hyperholomorphic) functions 

Min) = {f = fo + h6, + f262 eC\n,R'): Df = 0}. (1.2) 

The fact that M{n) C 'H(fi), where 'H(fi) = URiQ) x n^iQ) x n^iQ) is the set of M^-valued 
harmonic functions in Q, follows immediately from the factorization A = DD = DD of the 
Laplacian on M^. 

1.2 Antimonogenics and ambigenics 

We say that / is (left) antimonogenic when Df = 0. A monogenic constant is a function 
which is simultaneously monogenic and antimonogenic: Df = Df = 0, or equivalently, 
dof = df = 0. 

It is unavoidable that Df and Df need not take their values in even when / does. 
However, due to the fact that 

-esDf 63 = -63(^0 + <9)e3(-e3)(/o + /)e3 
= {do_-d){fo-f) 
= Df 

we have 

Df = O^Dj=0^fD = 0, 

and consequently. 

Proposition 1.1. A function is left monogenic if and only if it is right monogenic. The set 
of conjugates of monogenic functions 

M{n) = {J: feMm (1.3) 
coincides with the set of antimonogenic functions in Q. 

This is of course quite different from the situation for monogenic functions in H, where 
left- and right-monogenicity are different. The proposition allows us to write A4{Q) r\Ai{Q) 
for the set of monogenic constants. If / G fl then dof = df = 0. Thus a 

monogenic constant / does not depend on xq and can be expressed as 

/ = Co + /lei + /2e2 (1.4) 
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where Cq G M is constant and the quantity /i — 2/2 is an ordinary holomorphic function of 
the complex variable Xi + ix2- There are natural projections of Ai{Q) onto the subspaces 



ScM(fi) = {Sc/: feM{n)} c nu{n), 

VecM{n) = {Vec/: f e M{n)} C H{o}m<^). 

and by Proposition O we see that ScM{n) = ScM{n), VecM{n) = VecM{n). An 
element of + will be called an ambigenic function; its decomposition as a sum 

of a monogenic and an antimonogenic function is unique up to the addition of a monogenic 
constant. 

Consider -L2(fi, M'^) with the real-valued inner product (/, g) = Sc fg dV = J^{fogo + 
fi9i + f2g2)dV. Write 

M2{n) = M{n)nL2{n). 

It is somewhat inconvenient that A^2(^) is not orthogonal to A^2(^)- However, we have 
automatically that Sc A^2(^) -L Vec A^2(^) since any scalar function multiplied by cq is by 
definition orthogonal to any combination of Ci and 62- This fact gives us an orthogonal 
direct sum decomposition of the space of square-integrable ambigenic functions, 

M2{^) +M2{^) = ScA^2(fi) © VecA^2(fi)- (1-5) 

We will always assume that ^7 C is connected. 

Lemma 1.2. J77| / Suppose Q is simply connected, (a) Let fo G HrIQ). Then there exists 
f G A4{Q) such that Scf = fQ. This f is unique up to an additive monogenic constant. 

(b) Let /i,/2 G T-Lu{fl). A necessary and sufficient condition for there to exist f = 
/o + / G A4(f2) such that Scf = fo, / = /iCi + /2e2, is that 82 fi = dif2. When it exists, 
this f is unique up to an additive scalar constant. 

When the operator D is applied exclusively to scalar-valued harmonic functions, i.e. 

we see from (II. 5p and Lemma [1.21 that this operator splits naturally to give exact sequences 

^ M ^ HK(fi) %ScM^O, (L6) 

^keid ^nuin) AVecM ^0. (L7) 

Here kei d is two-dimensional over M, consisting only of polynomials h{x) = cq + ciXq. 

Note 1.3. One could equally well embed in HI differently, for example by considering 
{0} © C e, thus writing 

X = XiCi + X2e2 + XjCs = -e^x = ^261 - Xie2 + XqCs 
f = fiei + /2e2 + /scs = 63/ = -/2ei + /ie2 - /qCs 
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and using the operators 



D = eidi - 6282 - 6383, D = ei8i + 6282 + 6383. 

The relationship Df\x = Df\e^x imphes 

Df = ^ Df = 

<^=^ grad / = 0, curl / = 

and consequently this alternative embedding, which is used for example in [121 IS], is equiv- 
alent to the form we are using in this article. Our spaces correspond to 

the spaces of left- and right-monogenic functions / in the sense of In this context the 
equations defining monogenicity are also known as a Riesz system [71 [T71 [HI [22] . 



2 Homogeneous monogenics 

We will mostly work in the ball = < 1} C M^. The real vector space Hir(B^) fl 
L2(B^) is stratified into the subsets 'H^^\m^) comprised of real harmonic functions which 
are homogeneous of successive degrees n = 0, 1, . . . It is well known that the elements of 
'H^\m^) are polynomials of degree n and that these real linear subspaces are orthogonal 
with respect to the inner product (■,■). Further, every square integrable harmonic function 
on B^ has a unique expression as a series formed of elements of these sets. 



2.1 Basis for M2(M>^) 

Many schemes have been devised to construct bases for spaces of homogeneous monogenic 
functions of given degree, from the classical construction of Fueter to diverse applications 
of symmetric sums of products; see [5l [121 [IB] . For one way is to proceed as follows. A 
well known orthonormal basis of 'H^''(B^), n > 0, is the system of 2n -|- 1 solid spherical 
harmonics 

f/o",^r,...,[/„",i>i",...,v;", (2.1) 

where = r"'U'^, = r^V^, are defined in terms of spherical coordinates xq = rcos9, 
Xi = r sin 6 cos ip, X2 = r sin 6 sin ip via the relations 

U^ie^if) = Pn(cOS^), 

f/;^(^,^) = Pr(cos^)cos(m^), 

V::,{e,^) = Pr(cos^)sin(my.), m = l,2,---,n. (2.2) 

Here P„ is the Legendre polynomial of degree n and is the associated Legendre function 
is given by 
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We will need the following standard identities 

a-t')iPn+int) = (n + m+l)P-(t)-(n + l)tOW, (2.3) 
(l-^')'/'(^nTi)'W = PZ'i\t)-m{l-t')-'/hP-,{t), (2.4) 

(i-t'y^'Pn^iit) = i^{pz'2it)-pr\t)), (2.5) 

2mtPr+iW = (l-t')^/'(Pr+V(t) + (n + m + l)(n-m + 2)Pr+l'W),(2.6) 
(^_^ + l)p-^(t) = (2n + l)tP-(t)-(n + m)P™i(t). (2.7) 

The spherical harmonics (12.21) are polynomials when expressed in cartesian coordinates 
(xq, Xi, X2). One obtains a basis for the space A1(")(B3) of homo geneous monogenic functions 
of degree n formed by the 2n + 3 polynomials 

= Z^fc+'], m = 0,---,n + l 

= DlK^'l m = l,---,n + l, (2.8) 

These are monogenic by construction due to the factorization of the Laplacian. A detailed 
explanation of the analogous construction for — >^ EI is found in [S], [12]; the specific 
construction for — )■ M.^ given here appears in fTfi [18] and it is shown that X^, form 
an orthogonal basis for A^*^")(B^). For reference we recall that the derivation (12. 8 p produces 
the monogenic basis elements in the following form, 



= COS rriLp + (P^ COS Lf COS m^f — sin ip sin mi/)) ei 

= '^"^ sin "^V' + (-B^ COS (f sin my? + sin cos my?) Ci 

+ (P^ sin sin mcp — cos ip sin m(p)e2 ] , (2.9) 



t=COS ) 



where 

= i((l-t^)(Pr+i)'(t) + (n + l)tP™i(t))|*=cos., 

= ^(v^r^i(^r+i)'w-(^ + i)v/r^^r+iw) 

777 

C^;^ = ^^7^^n"ll(^)k=cos.. (2.10) 

(We observe that other authors have used notation such as X^'^, y^'"!" where we write X^, 
Y^, in order to stress that these functions are defined in rather than on S"^.) Finally in 
L2(B^,M^) we consider the norm II/II2 = a/ (/, /), where the scalar product is over B^ as 
in section [L2] (rather than over S"^ as is the case of some other authors). The norms of the 
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solid spherical harmonics and orthogonal monogenic functions are given by 

mh = 



An 



(2n+ l)(2n + 3) 



\Ulh = IIK:II2 = J ,, ^Jf^^'^t (2.11) 

[Zn + lj{2n + dj [n — mjl 



, II2 



7i{n + 1] 
2n + 3 



ivr^W - - / 7r(n + l)(n + m + l)! 

\^mh - \\yj\2 - \l2{2n + 3){n-m + l)V ^^'^^^ 



when m > 1. 



The following explicit representation of the basis elements of A^'^")(B^) in terms of spheri- 
cal harmonics is stated in [18j without proof. Since our results depend on this representation, 
we will prove it here in detail. 

Theorem 2.1. Write 

^ (n + m)(n + m + l) 
Cm = ^ • 

For each degree n > 1, the basis elements for the homogeneous monogenic polynomials of 
degree n are given by 

Y'^ — (""- + ""^ + 1) j}n _ ( ^rijjn _ ^ ffn ) „ i ( \}n i T>« \ p 



~ 2 ('^m^rn-l j ^1 + ) ^2, (2-13) 



where 1 < m < n + 1 . 

Proof. Recalling f ll.6p .f lL7|) we see that since and V^'^^ are scalar valued, the definition 
(12. 8 p may be expressed as 



so the components of X^, given in (12. 9p are precisely the partial derivatives of W!^^ , Kn"*^^- 
Consider the formula proposed for Xq in (12.130 : i.e., assume for the moment m = 0. By 
(12. 9p . it is necessary to prove the three equalities 

^^f/o" = A^, (2.14) 

hj^ = B^cos^, (2.15) 

^V,^ = B^sin^. (2.16) 



Substitute the definitions (12 ■2p and fl2.10p into these equations. We see that (12.141) is given 
immediately by (12. 3p . To verify equations (I2.15p . (I2.16P we must prove that 

(1 - t2)i/2^(P„+i)' -{n + l){l- ef^Pn+i - Pi = 0. 

where for brevity we write Pn+i in place of P„+i(t). The left hand side may be expressed as 

(1 - t'Y/\t{p^+,y -{n + i)p„+i - (p„)') 

and when we substitute the values for (Pn)', (Pn+i)' given by (12. 3p we the result is zero 
according to Thus (12151) and (I2A6|) hold. 

Similarly the desired formula for X^, l<m<n+l, is equivalent to proving the three 
equalities 

Al^cosmip, (2.17) 
B'^ cos ip cos rrup — Cl^ simp sin mip, (2.18) 
P^sin(y9cosm(y9 + C^cosv^sinmy). (2-19) 

As before, substituting the definitions (12. 2p and (I2.10p we obtain again that (I2.17P reduces 
immediately to (12. 3p . For equation (I2.18P we equate the coefficients of cos ip cos rmp and of 
siny9sinm(y9 on both sides; for (I2.19P we use the coeficients of sinipcosmip and cos(p sinrrvp. 
From this it is seen that proving equations (I2.18P and (I2.19P is equivalent to proving 

+ m)(n + m + l)P™-i - Pf^^ + 2(1 - t^)^/^t(P™ i)' 

-2(n + l)(l-t2)i/2p-^ = 0, (2.20) 
(n + m)(n + m + l)P„"^-^ + Pf+i-2m(l-t2)-i/2pn^i = 0. (2.21) 

Let us verify (I2.20p . Eliminate the derivative (P^^)' via (12. 3p to convert the left hand side 
into 

{n + m){n + m + l)Pf - P™+i + 2{n + m + l)t{l - t2)-i/2pm 
-2(n + l)(l-t^)V2p-,. 

We note that (12. 6p provides a value for {n + m){n — m + 1)P^~^, so we are led to decompose 
n + m + 1 = 2(?2 + l) — (n — m + 1) and to arrange the terms as follows, 

2{n + m)(n + l)P„"^-i - (P„"^+i + {n + m){n - m + l)Pr"') 

+ 2(1 - f )-^/\in + m+ l)tP^ - (n + l)Pr+i), 
= 2{n + m)(n + l)P™-i - (1 - e)-^'\{n - m + l)tP™ - 2{n + 1)P„™ J. 

Multiply this by (1 — t^)^/^, and then substitute in the first term the value for (1 — t^)^/^P^"-^ 
provided by (12. 5p to arrive at 

{n + m)(P- 1 - Pr_i) + {2n + l)tPr - {2n + 1)P™ ^ 



{n + m + 1) 



Jjn I Tjn 
m m—1 ' ^ m+1 
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which clearly vanishes by { \2.7\i . thus proving fl2.20p . 

Now we verify equation (I2.2ip . Its left hand side can be written as 

2m{n + m)P;;'-' - 2m(l - tY'^^P^+i + + {n + m){n-m + l)Pr~' 

and by applying fl2.6p to the last two terms it is equal to 

2m(n + m)P™-^ - 2m(l - t2)-l/2(prn ^ _ ^prny 

We may divide by 2m(l — t^y/'^. Using fl2.5p (with — 1, m — 1 in place of m, n) to replace 
the first term, the result is zero by (12. 7p . This proves fl2.2ip and thus the formula for X^. 
The formula for hkewise reduces to (I2.20p and fl2.2ip . so (12.130 is verified for all cases. 
□ 

Theorem 12. II immediately yields a basis for Vec A^, which turns out to be orthogonal: 

Proposition 2.2. For each n>{), the set 

{ VecX"^ = 9f7;^+\ < m < n + 1} U { FecyjJ = at>+\ 1 < m < n + 1} 

is an orthogonal basis for yecA^''"^(B^). The union of these sets over all n > is an 
orthogonal basis for VecA4(M^) , and the norms of the basis elements are given by 



VecX^l 



TTn{n + 1) 



(2n + l)(2n + 3) 



x-nii II A^nii2 / vr(n2 + m2 + n)(n + m + 1)! 

VecXZW = WVecYJ^r = W 7777^—7^—^^7 TTTT ' - ^• 

2{2n + lj{2n + Sj{n — m + 1)1 



' m I 



Proof. The given set is clearly a basis of Vec A^^"^(B'^) because of (11.71) . We need to see that 
it is orthogonal. Choose elements f,g in the basis for the monogenics {Xq, X^, Y^}, and 
express them as f = fo + f , g = go + g, with scalar parts /o, go- Clearly 

{f,9) = {fo,9o) + {f,9)- 



By Theorem 12. H the scalar parts of the monogenics run through (scalar multiples of) the 
spherical harmonics, and thus are orthogonal. Suppose f ^ g. Then (/, (?) = and (/o, go) = 
0, so (/, ^) = as desired. To calculate the norms, now suppose f = g. Then (/, /) = 
12 - ||/o||2. Thus if / = Xg", then by (Km . fl212|) . and (12^3]) 

7r(n + l) {n + iy An _ TTn{n + 1) 

\J 1 J) 



{2n + 3) 4 (2n + l)(2n + 3) (2n + l)(2n + 3) ' 

Similarly, if / = X^ or / = 1^ for m > 1, 

7r(n + l)(n + m + 1)! (n + m + 1)^ 2n{n + m)\ 



ifJ) 



2{2n + 3){n - m + 1)\ 4 (2n + l)(2n + 3)(n - m)! 

7r(n^ + + n)(n + m + 1)! g 



2(2n + l){2n + 3)(n - m + 1)! ' 
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We now give a basis for the ambigenic functions. It must be noted that the monogenic 
constants X'!^^-^ and Y^_^_-^ are the negatives of their conjugates, so care must be taken to 
count them only once. Thus the following functions thus form a (not orthogonal) basis for 
the ambigenic functions on B^: 

vn vn x^n x^n \rn -yn y-n vn Vn Vn Vn vn (r\ r\r\\ 

y\Q , . . . , y\.„, y^n+l^ 1 5 • • • 5 -' ra 5 ^0 ' ^1 ' ' ' • ' ' 1 ' • • • ' ra ' -'n+l' 



Proposition 2.3. Let n > 0. The 4n + 4 functions 



m 

m 
vn,- 
m 



m = 0, 



n + l, 



m = 1, . . . , n. 



X" — n'^ X'^ m = Q 



m = 1, . . . , n + 1, 



(2.23) 



where 



n — 2m + 1 



(0<m<n), <+i = 



- (n + l)(2n + l) 

/orm an orthogonal basis for the space of square integrable ambigenic functions on which 
are homogeneous of degree n. 

Proof. Since these are 4n + 4 ambigenic functions, it suffices to prove the orthogonality. 
First we calculate the scalar products of each X^, and its conjugate. By Theorem 12.11 
and (12. lip we obtain 



n{n + l) 



(2n + l)(2n + 3)' 

TT{n - 2m^ + l)(n + m + 1)! 



2(2n + l)(2n + 3)(n-m + 1)!' 



(2.24) 



for 1 < m < n + 1. Since the set {Xq, X^, Y^: m = 1, . . . , n + 1} is an orthogonal basis of 
A^(")(B3), it follows at once that 

-,^1 I - -,^1 )-\Jm J^l ) — ■> I / " U- 

Further, by f l232|) 

if m = / = 0, 
if m = / ^ 0, 



' 7r(n + 1) 

2/2 + 3 

7r(n + 1)(?7, + m + 1)! 
2(2n + 3)(n-m + 1)! 
ifm^^/. 



Substituting fl2A2ll and ( 12:241) we find that 
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for all l,m, and also that 



4:7m{n + 1)^ 



(2n + 3)(2n + 1)2 

27r(?7,2 + + ?T,)(?T, + m + l)(n + m + 1)! . 



ii m = I = 0, 



{n + l){2n + 3)(2n + iy{n - m)\ 







if m = / ^ 0, 
if m 7^ Z. 



The calculation of the scalar products for {1^}, as well as for the mixed cases, is similar. 
□ 



3 Contragenic functions in 

3.1 Existence of contragenics 

From now on we will abbreviate M^") = A^(")(B3), AT^"^ = M^'^\m^), = 'H^^\m^). 
Because of the correspondence of monogenic constants with holomorphic functions described 
in section [H and since the real homogeneous harmonic polynomials of degree n in the complex 
variable x + iy are linear combinations of Re {x + iy)^ and Im (x+iy)"', the space 
is 2-dimensional over M for n > 1. We summarize in Table 1 the dimensions over M of the 
relevant spaces of functions in B^. Recall that T-L denotes R'^- valued functions. 



Space of polynomials 


n = 


n > 1 


-T/(") 


1 


2n + l 




3 


2n + 3 




3 


2 




3 


An + 4 




3 


6n + 3 



Table 1: Dimensions over M of spaces related to monogenic polynomials in B^. 

From the last two rows of this table comes the following notable fact. 

Theorem 3.1. Not all L2-harmonic functions are ambigenic: M.2+M.2 is a proper subspace 
ofn{M^)nL2{M^). 

Definition 3.2. In any domain Q, a harmonic function h G HlQ) fl L2{ft) is called con- 
tragenic when it is orthogonal to all square- integrable ambigenic functions, i.e., if it lies 
in 

Afin) = {M2in) + M2in))^ 



11 



where the orthogonal complement is taken in T-L{Q) fl ^2(^2). In we have the orthogonal 
complements in the spaces of homogeneous harmonic polynomials 'H^^\ 

As observed in the introduction, there is no direct analogy to be found for contragenic 
functions in hypercomplex analysis on C or EI or more general Clifford algebras since in 
those contexts all harmonic functions are known to be ambigenic. 

Let n > 1. It follows from that A^(")(fi) nM^''\n) C VecA^(")(fi). By Lemma 
ll.2[ when Q is simply connected we have Sc A^*^"'*(f2) = 'H^\Q). Returning to B'^, we have 
specifically dimScAl*^") = dim'H^-' = 2n + 1. From (11 .Sp it follows that dimVec A^*^") = 
2n + 3. (Since this is equal dim A^*^"\ this means that given the vectorial part /, the scalar 
part /o is uniquely determined.) Thus 

dimAT^") = 2n - 1 (3.1) 

when n > 1, while dimA/''-^^ = 0. 

The following result is a simple consequence. 

Theorem 3.3. Let h E M = A/'(B^). Then h can he uniquely expressed as a sum in L2 

00 

n=l 

where /i^") G TV^") . 

Note 3.4. It is easily checked by a dimension count that the analogue of Theorem 13.11 for 
"clasical" monogenic functions H — t- EI does not hold; i.e., all harmonics are ambigenic over 
H. For n > one has that the homogeneous monogenics in B^ form a right vector space 
over EI of dimension ^{n + l){n + 2) |2B], and the same is true for the antimonogenics. The 
monogenic constants have dimension n+1 over EI (See example [2j). Since the dimension of 
the harmonics from EI to EI is (n + 1)^ over EI [28], it follows that every harmonic function 
can be expressed as a sum of a monogenic function and an antimonogenic function. 

It may also shed light on the situation to see what fails when one attempts to express a 
harmonic M^-valued function in terms of monogenics and antimonogenics. For scalar-valued 
/o G 'Hr(^^) there is in fact no problem, since 

/o = ^(/o + ^0) + ^(/o - ^0) 

where /o + ^0 is the completion of /o to a monogenic function as given by Lemma 11.21 
When we are given a general harmonic function /q + /iCi + 7262 G 7^(1^) and complete each 
component separately, we obtain analogously 

/o + fiei + /2e2 = ^((/o + ^0) + (/i + gi)ei + (/2 + ^2)62) 

+^((/o - ^0) + (/i - 9i)ei + (/2 - ^2)62). 
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However, the two functions added on the right hand side of this equation, while in the kernels 
of D and D respectively, need not take their values in M^. Thus we do not obtain in this 
natural way a representation as an ambigenic function. 



3.2 Basic properties of contragenics 

Although our main interest is in B'^, we observe also some relations which hold in more 
general domains. Consider an arbitrary contragenic function h = /loCo + Ziiei + /i2e2 € 
By the orthogonal decomposition (II. 5p . h ± Sc A^2(^); which implies {ho, ho) = so in fact 
Hq = 0. Thus h is of the form 

h{xo, xi, X2) = hi{xo, xi, X2)ei + /i2(a;o, Xi, X2)e2. (3.2) 

In particular, contragenic functions are never invertible. Equation f ll.Sp also gives the prop- 
erty h _L VecA^ 2(^)5 which by Lemma [1.21 can be expressed as 

/ {hifi + h2f2)dV = (3.3) 
Jn 

whenever ^1/2 = 82/1 with fi and /2 harmonic. 

Work on monogenic functions in different contexts has focused on standard domains such 
as spheres, ellipsoids, cylinders, rectangles, etc. (see for example [E]). Let Qi C be a 
domain which enjoys the symmetry that {xq,X2,Xi) G Qi whenever x = {xo,Xi,X2) € Qi. 
For f : Qi -^M write 

f*ix) = fo{Xo,X2,Xi) + f2{Xo,X2,Xi)ei + fl{Xo,X2,Xi)e2 - f3{Xo,X2,Xi)e3 

for X E Qi. Thus (/*)* = /, and it is easily seen that 

if 9)* = 9*r (3.4) 

for f,g:Qi^ M^, and 

Din = (DfT (3.5) 
when / is different iable. From this we have 

Proposition 3.5. The involution f ^ f* preserves and J\A{Qi). If is simply 

connected, the involution preserves as well. 

Proof. Let / G Mini). By ([33]), D{f*) = 0* = 0, so /* G A^(fii). Similarly M is invariant. 
Now let h e A/'(^^l). For g e HM(f]i) n ^2(1^1) take / = dig ei + d2g 62 G Vec A^2(^^i)- We 
find that 



{h*, f) = I {h2{xo,X2,xi)dig{x) + hi{xo,X2,xi)d2g{x)) dV 

{h2{xo,X2,Xi)d2g*{xo,X2,Xi) + hi{xo,X2,Xi)dig*{xo,X2,Xi)) dV = 
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by (I3.3p . where g*{x) = g{xo,X2,xi). Assuming fli is simply connected, it follows from 
Lemma [1.21 that / ranges over all of Vec A^2(^i); thus h* G Af{^i) as claimed. D 

We conjecture that the simple connectedness is not necessary for the invariance of A/'(fii) 
in Proposition 13. 5[ 

In the unit ball we have the following characterization of contragenic functions via inte- 
gration over the unit sphere. 

Proposition 3.6. Let h G H^'^^M^), h = hei + /laCs. Then h G A^(")(B3) z/ and only if the 
equality 

I hig dxo A dx2 = / h2gdxoAdxi (3.6) 

holds for all g G ni^\M^). 

Proof. We have the following relation of differential forms, 
{hidig + /i2<925') dxQ dxi dx2 = 

d{—hig dxo A dx2 + h2g dxQ A dxi) — {dihi + ^2/12)5' dxo dxi dx2. 
Since harmonic functions of differing degrees are orthogonal over B^, integration leaves 

{h, dg) = / d{—hig dxo A dx2 + h2g dxo A dxi) 

which is equal to whenever h G A/", and Stokes' theorem gives (13. 6p . Conversely, if (13.61) 
holds, then {h,dg) = 0. Since dg ranges over all of VecAi^^\ h is orthogonal to A^*^"^. 
Further, since h is trivially orthogonal to A/l^'"^ for m 7^ n, it is contragenic. D 

Corollary 3.7. Let h G H{W^), h = hid + h2e2. Then h G 7V(l3) ^^ly ^) /^^/^^ 

for all g G Hr{W^). 

Proposition 13.61 and Corollary 13.71 together with the bases given in section Ej provide 
an algorithmic method for determining when a given harmonic function, expressed as a 
convergent series in HI fl L2, is contragenic. However, it is not likely that there is a simple 
characterization of contragenics purely in terms of derivatives for general Q, and indeed we 
know of none even for B'^. 



3.3 Bergman kernel for \ecA4 

The natural generalization of the holomorphic Bergman kernel [H] from the context of 
holomorphic functions in C to that of monogenic functions in EI is described in [23 ES]- A 
generalization for functions in was defined and studied more recently in [9]. We restate 
some of the main facts in the present terminology, and then give a new "Bergman kernel" 
which is more appropriate to the subject at hand, as it provides another characterization of 
contragenic functions. 

The following result establishes that evaluation at a fixed point is a continuous linear 
functional on VecA4(f2); it suffices to work in B^ for the basic estimate. 
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Proposition 3.8. Let f e VecM. Then for all x G M^, 



\fm<c\\fh 

where C = y/3/{An). 

Proof. Since / is liarmonic, using tlie ortliogonal basis of Proposition 12.21 we see that the 
constant /(O) is orthogonal to / — /(O) in L2{M^). (More simply, one may just observe 
that the constants are orthogonal to all harmonic functions which vanish at the origin, a 
statement of the Mean Value Property of harmonic functions.) Therefore 

\fm' = c'[ \m\'dv = c\\\f\\i-\\f-fmi)<c'\\f\\i □ 



The underlying idea, given a closed subspace L2 and an orthonormal basis {^Pk} of 
A, is to form an integral kernel B{x,y) = J2kVk{x)Lpk{y) which automatically enjoys the 
reproducing property f{x) = J B{x, y) f{y) dy for f E A, and projects L2 orthogonally onto 
A. However, given an orthornormal basis {(pk} of A^2(^); there is a problem if we try to 
construct a Bergman kernel in this way because MI2 is not closed under multiplication: the 
integrand of 

{Bo{x,y)eo + Bi{x,y)ei + B2{x,y)e2 + B3{x,y)e3) ■ (/o(y)eo + /i(?/)ei + /2(l/)e2) dVy 

contains the term (-B1/2 — -B2/1 + -83/0)63, whereas / should be M^-valued. Thus one needs 
an additional condition 

J (i?i/2 - 52/1 + 53/0) dy = 0. (3.7) 

In [9] this is dealt with by working in the subspace of L2 corresponding to functions / for 
which this property holds. We will take a different approach here, constructing a Bergman 
kernel for Vec A^2 rather than A^2; and not requiring a special condition such as 03.71) . 

The scalar product restricted to L2{Q, {0} © © {0}), a Hilbert space which contains 
VecA^2(^^) closed subspace, is 



(/, 9) = J if 191 + f292) = -Sc J fg 



since / = — /. Take an orthonormal basis {tpk}kLi of Vec A^2(^) over M, and write ipk = 
i'k,!^! + 'ipk,2e2- Define the following functions Q x Q ^ {0} © © {0}, 

00 

bn,i{^^y) = -^^k,i{^)^k{y), 

k=l 

00 

bn,2ix,y) = -^^k,2ix)^pk{y)■, (3.8) 



k=l 
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it can shown by means of Proposition 13.81 that these series converge uniformly on compact 
subsets of r2 X r2. 

Definition 3.9. The Bergman operator for VecAi{fl) is defined by 

Bn[f]ix) = Sc ( / bn^i{x,y)f{y)dVy) Ci + Sc ( / bn^2{x,y)f{y)dVy)e2 



for all / e L2{n, {0} © ^ {0}) and x e Q. 

It is shown in the traditional way that bQ^i{x,y) and bQ2{x,y) are independent of the 
orthonormal basis chosen. Since we can express elements of Vec f = J2 (^ki^k {cik G K), 
the following reproducing property is easily checked. 

Theorem 3.10. The linear operator Bq projects L2{^1, {0} © © {0}) orthogonally onto 
VecM2i^)- particular, Bnlf] = f for f e VecM2i^). 

Note that two separate integral kernels (13.81) are necessary in the definition of the operator 
Bq because the scalar product on Vec Ai (Q) is only bilinear over the reals. 

For B^, in terms of the specific basis := dU^, := dV^ of VecA^*^") given in 
Proposition 12.21 one can define kernels for each degree n, 

n n 

bUx,y) = -^u:;^i(x)u:;(y) - ^<^i(a;)t;;;(y), 

m=0 m=l 
n n 



m=0 m=l 



and then can form operators S^"-* analogously to Definition 13. 9[ These operators project 
the harmonic functions onto VecA^*-"\ and the Bergman operator B for B'^ is their sum 
B = Yl'^=o B^'^^ ■ It would be interesting to express ( 13. 9p in closed form. 

The following is an immediate consequence of the foregoing, and with the formulas (13.91) 
allows one to detect computationally when a harmonic function is contragenic, or close to 
contragenic in the L2-sense. 

Corollary 3.11. Let h = hiCi + /i2e2 G ^2(^2) he harmonic. Then h G N'{VL) if and only if 



4 Construction of homogeneous contragenic polynomi- 
als 

In this section we will give an explicit construction of a basis of A/"^"-* for every n = 0, 1, . . . 

One possible approach would be as follows. The basis for M.^'^'' + Ai^^^ given in Propo- 
sition 12.31 is orthogonal, so one may extend this to a basis of "H^"^ by choosing suitable 
linearly independent triples of spherical harmonics and applying the Gram-Schmidt process 
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to produce the contragenic polynomials of degree n. However, this procedure is quite costly 
numericallyj§ and leads to little insight regarding contragenic functions. 

Here we give a direct construction of the contragenic homogeneous functions. From Table 
lOit is clear that A^^^) = {0}. 

Theorem 4.1. Let n > 1. Write d^ = {n — m){n — m + l). The 2n — 1 functions 

:= V,^e,-Ul'e2, 

: + Ul+i)ei + + C+i)e2, (4.1) 

for 1 < m < n — 1, form an orthogonal basis for Af^"'^ over M. 



Proof. First we show that the functions (14. ip are contragenic: it is sufficient to show that 

each one is orthogonal to '•"^ U " . As we have already noted, since they have no scalar 
parts it suffices to show that each one is orthogonal to VecA^^'^^ and to do this, by (11.71) 
we may use the basis VecA^'-"^ obtained by dropping the scalar parts of the basis for A^^"-* 
given by Theorem 12.11 Let 1 < m < n — 1 . From Theorem 12.11 

(Z;;,^,VecXo") = l{{dlK_,,U]:) + U-) + {dl^^^ 

By the orthogonality of the spherical harmonics and V^, this scalar product is equal to 
zero. Next we observe that for 1 < k < n, 

+ {d^U^_i — '^V^k-i + ^^fc+i) 

= 0, 

since U'^^^ is orthogonal to V^^^. 
Finally, it remains to check that 

= 0. 



^Calculations in Mathematica on a desktop computer with 4 Gb of RAM have saturated the memory 
when attempting to calculate contragenic homogeneous polynomials of degree n > 7 via Gram-Schmidt as 
described here. 
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Once again it is immediate that this is true under the condition k m, k m + 2 and 
k ^ 171 — 2 since all of the scalar products involved vanish. Now consider k = m. Substituting 
the equation (12.111) we obtain that 



{n + m){n + m + 1) 27r(n + m — 1)! 



—2(n — m)(n — m + 1) 
4 



4(2n + l)(2n + 3) (n-m + 1)! 
2 -271 (n + m + 1)! 



4(2n + l)(2n + 3) (n - m - 1)! 
0. 



When = m + 2, we see by (12. lip that 

= 

and the case = m — 2 is similar. Therefore {Z^ ^, Vec Yj^) = for all k. This completes the 
proof that ^ is contragenic. The proofs that Zq and _ are contragenic are analogous. 

Now we show that these functions form an orthogonal basis. Since dim^TV^") = 2n- 1, 
it suffices to show that they form an orthogonal collection. The only nontrivial cases are 
(Zq,Z£^), (Z^ Z^_|_2^+), (Z^^+, Z^_2,+), (Z^_,Z^+2,-) (Z^_,Z^_2_), all of which 
are seen to be zero by repeated applications of (12. lip . D 

Corollary 4.2. The set 

{Zl, Zl^^, l<m<n-l, n>l] 
is an orthogonal basis for M = A/'(B'^). The norms of the basis elements are 



I 

1^0 I 



I 7^ I 



A7m{n + 1) 
(2r2 + l)(2n + 3)^ 



^7r(n2 + + n){n + m — 1)! 
(2n + l)(2n + 3)(n-m-l)!' 



/or 1 < m < n — 1 . 

Proof. It is only necessary to establish the values of the norms. By (12. lip 

7n\\2 _ I ryn ryn\ _ /f>n f>n\ , /f>n f>n\ _ 47rn(n + 1) 



\zn' = {Z-, ZD = {v,\ vn + (t/r, un 



(2r2 + l)(2n + 3) 
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and, 

II yn ||2 / 7" 7*^ \ 

ll^m,+ ll ~ \^m,+5 ^m,+/ 

47r(n — rnf'{n — m + l)^(n + m — 1)! 47r(r;, + m + 1)! 

(2n + l)(2n + 3)(n-m + l)! ^ (2ra + l)(2ra + 3)(ra - m - 1)! 
87r(n^ + + + m + 1)! 
(2?T, + l)(2n + 3)(ra + m)(n + m + l){n — m — 1)\ 
'&'K{n^ + + + m — 1)! 
(2n + l)(2n + 3)(n-m-l)!' 

The calculation for _ is similar. D 

Note 4.3. The involution / h-t- /* of Proposition 13.51 sends ^ to (a multiple of) _ for 
some, but not all m. 

5 Conclusions 

Consider a triple / = /oeo + /iei + /2e2 of harmonic functions in a domain VL. We have shown 
that / has a natural decomposition f = g + h where g is ambigenic and h is orthogonal in 
L2(f2) to all ambigenic functions. The existence of nontrivial contragenic functions raises 
the following question. Suppose that Q has smooth boundary and / is defined only on 
dfl. When is the harmonic extension of / to the interior of Q monogenic, ambigenic, or 
contragenic? How do the boundary values of the monogenic, ambigenic, or contragenic part 
of the extension relate to the original /? 

Further, it remains to investigate bases of contragenic functions in domains of other 
than B^, as well as analogous notions of contragenicity with respect to other scalar products, 
for example in weighted inner product spaces or with respect to the Fischer product [23t l27]. 
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